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We investigate doping of a two-orbital chain with mobile S = 1/2 fermions as a valid model 
for Y2- x Ca x BaNi05. The S = 1 spins are stabilized by strong, ferromagnetic (fm) Hund's rule 
couplings. We calculate correlation functions and thermodynamic quantities by DMRG methods 
and find a new hierarchy of energy scales in the spin sector upon doping. Gapless spin-excitations are 
generated at a lower energy scale by interactions among itinerant polarons created by each hole and 
coexist with the larger scale of the gapful spin-liquid background of the S = 1 chain accompanied 
by a finite string-order parameter. 



Since Haldane's famous conjecture about the existence 
of a spin-gap in one-dimensinal integer Heisenberg (HB) 
spin-chains the intriguing physics of S = 1 chains 
and 5=1/2 ladders has been investigated in numer- 
ous theoretical and experimental studies and especially 
hole doping of a spin-liquid is of great importance for 
the Mott-transition in the context of high- Tq supercon- 
ductivity. However, doping of a Haldane 5=1 chain 
with mobile holes has only recently been realized in the 
material Y2BaNiOs 0. In the undoped case, the two 
active Ni 2+ orbitals are coupled by a strong fm Hund's 
rule coupling Jh, giving an almost ideal 5 = 1 HB chain 
with a large spin gap of Ag = i w lOOif and very weak 
inter-chain interactions. By substituting Ca 2+ by Y 3+ , 
mobile holes can be introduced in the chain. The most 
remarkable experimental features are a reduction of the 
resistivity pd c by several orders of magnitude, and the 
creation of states with 5 between 1 and 3/2 per impu- 
rity inside the Haldane gap Q . Further the temperature 
dependence of the resistivity no longer follows thermal 
activation over a barrier. 

While doping of 5 = 1 HB chains by static as well as 
bond impurities UH has been studied by several authors, 
only little is known about doping of mobile 5 = 1/2 
fermions in 5 = 1 chains. A first theoretical investiga- 
tion limited to the case of weak hopping amplitudes || 
found states inside the Haldane gap with 5 larger than 
5=1/2 per impurity, whereas other studies based on 
the weak coupling approach find that the spin gap is im- 
mediately destroyed upon doping for models with a level 
difference and fully mobile electrons . Numerically, the 
spin dynamics of mobile holes in Y2BaNi05 have been in- 
vestigated by exact diagonalization of small finite chains 
for effective one ||] and two band models Jl|], and based 
on the assumption of infinitely strong Hund's rule cou- 
plings Jh states below the Haldane gap have been found. 

Among the most interesting properties of these sys- 
tems is a competition between fm ordering induced by 
the double exchange (DE) mechanism and antiferromag- 
netic (af) order due to direct exchange. Other problems 
are the persistence of the spin-gap and the characteri- 
zation of the doped phase. In this Letter we would like 
to address for the first time this competition in an un- 
perturbative manner, with finite values of the Hund's 
rule coupling and fully mobile electrons for system sizes 



significantly exceeding the correlation length of the un- 
doped system. We do so by investigating groundstate 
properties and various correlation functions by the den- 
sity matrix renormalization group (DMRG) method [ jTo| , 
and thermodynami c q uantities by the thermodynamic 
DMRG (TDMRG) fjlf. This reveals a new hierarchy of 
energy scales in the spin sector with gapless spin-modes 
given by the polaron-polaron interactions and a second, 
much higher lying energy scale stemming from the gap- 
ful spin-liquid background of the 5 — 1 chain. We also 
give first numerical evidence that the spin-gap is imme- 
diately destroyed upon doping of mobile electrons in a 
5=1 chain with level difference. The hierarchy struc- 
ture found for doping of static holes into the spin-gap 
background |l2| ] surprisingly persists even in the metallic 
phase with fully mobile hole doping. 

We investigate a model with strong Hund's rule cou- 
pling Jh < and a fully occupied lower band, as the 
d x 2_ y 2 orbibtal of Ni 2+ in Y2_ x Ca x BaNi05 is almost lo- 
calized. The results for a model with equal doping in 
both bands will be published later Further we as- 

sume af couplings J > between neighboring sites in the 
upper (3d z 2_ r 2) orbital and Jd > between neighboring 
sites in different orbitals. Thus our Hamiltonian reads 



J ^2 (Sj,iSj + x,i - -n jt xnj + i t i 

3 x 



where the first index of the double index j, i denotes the 
lattice site and the second index the orbital, the projec- 
tion operator V prohibits double occupancy of a site in 
order to take account of the strong on-site Coulomb re- 
pulsion and the rest of the notation is standard. If not 
otherwise mentioned we always set J = Jd and assume 
J H > J,Jd,t, typically -J B = KM = 20 J = 20 J d - We 
have also performed calculations with two different chem- 
ical potentials for the upper and lower band, but we find 
finite hole densities in the lower band only at the high- 
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est temperatures T > t in that case, and the results are 
identical to the above case for T < t. 

The TDMRG allows us to calculate thermodynam- 
ical properties in the thermodynamic limit of infinite 
system size. In general we use finite Trotter-time step 
sizes of At = Q.2t and extrapolate from the grand 
canonical ensemble with fixed chemical potential to con- 
stant particle density first. We keep between 60 and 80 
states per system and environment block, and use a re- 
biorthogonalization method fl4|| if we encounter numer- 
ical instabilities. With the groundstate DMRG we have 
considered system sizes of up to 2 x 256 sites and kept 
up to 1300 states per system and environment block. 

We begin by considering two limits which are easy to 
understand. The first case is J = Jd = 0. Away from half 
filling, all spins are ferromagnetically aligned by the DE 
mechanism Jl5| in order to gain kinetic energy. The sec- 
ond case is at half filling, where upon switching on the af 
couplings J, Jd, the system can be mapped to an af 5 = 1 
HB chain, with a spin-gap and exponentially decaying af 
spin-correlations. For any finite value of J, Jd, and finite 
hole doping we will therefore have competition between 
fm order induced by the DE mechanism and the spin liq- 
uid state of the 5 = 1 chain. Let us briefly sketch the 
picture that emerges from our numerical calculations, be- 
fore investigating the physical quantities in detail. Each 
hole doped into the 5 = 1 chain is surrounded by a small 
fm cloud created by the DE mechanism, and we will call 
this a polaron in the following. However, in the parame- 
ter regions we studied in this paper, the polaron is only 
a weak local perturbation of the spin- liquid background, 
and the correlation length remains close to its original 
value. Among the polarons, Akp and 2\zf charge den- 
sity wave (CDW) order is then stabilized, and the low- 
est lying spin excitations are weak magnetic interactions 
among the polarons. Due to these interactions, there is a 
spin-singlet groundstate already for very weak couplings 
J, J d > 0-2t and J H = -lOt. 

Let us turn to the detailed investigation of physical 
quantities next and address the question of whether the 
spin-gap is robust upon doping of holes first. The map- 
ping to the af 5 = 1 HB chain at half filling gives effective 
couplings J e = 3/4 J of the original 5=1 chain by per- 
turbation theory, hence we expect a spin gap of A s = 
3/4As =1 of the original value Ag =1 « 0.41J. In good 
agreement we obtain A s w 0.153* ± O .OOSt = 0.3 06J ± 
0.01 J from the TDMRG by lim T _o y/2Tc v /(3x) -> A s 
fusf . Alternatively, we have calculated the spin gap 
by finite size scaling A s — \miT^oA s (L;N — Ln), 
where A S (L;N = Ln) = A S (L;N) = E {L;N;S Z = 
1) - E {L;N;S Z = 0) and E (L;N;S Z ) is the ground- 
state energy of the system with N particles on L sites 
and total S z component of the spin along the z-direction. 
However care needs to be taken in order to avoid free 
5 = 1/2 spins at the end of the open chain. To keep 
the free spins out, we have imposed the boundary condi- 
tion that the upper orbital on the boundary sites remains 
empty. The result for different values of the hole con- 
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FIG. 1. Magnetic susceptibility \ for 

— Jh = 10t = 20 J = 20 Jd- Inset: Effective Curie-constant 
S(S +l)/3 per hole. 



centration in the conduction band is listed in Tab. |l| 
and agrees well with the TDMRG. Here we see that the 
spin-gap is immediately destroyed upon doping, and by 
noting that the smallest momentum in the finite open 
chain is fc m j n = n/L we can calculate the spin-velocity 
v a from A S (L) = w CT A: m i n . Thus we can determine the 
T = susceptibility by x° = where K a = 1 

due to the SU(2) symmetry. Both results are listed in 
Tab. |. We compare these findings with the TDMRG 
results for x{T) in Fig. |l|. The strong enhancement of 
x(T) at low temperatures suggests the formation of larger 
magnetic moments upon doping, in agreement with the 
T = results. In the inset of Fig. [j] we estimate the 
effective Curie-constant 7 = 5(5 -I- l)/3 created by each 
hole by subtracting the background of the 5=1 chain 
7 = (x — (1 — n'f l )xs=i)T /n 1 ^. The maximum value of 
7 « 0.5 near T = 0.3t is nearly independent of the hole 
doping and corresponds to 5 c ff ~ 0.7, in rough agreement 
with experimental findings H17L By further lowering 
the temperature, af interactions among the polarons re- 
duce 7 and finally lead to a singlet groundstate at T = 0. 
We have checked the rotational invariance by the DMRG 
from = S Z S Z — \/2{SfSj) which vanishes for a 

singlet groundstate and confirm this result within the nu- 
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TABLE I. The spin-gap A s for different values of hole dop- 
ing n c h in the conduction band and — Jh = lOt = 20J = 20Jtj. 
Also listed is the spin-velocity v a and the susceptibility %° a t 
T = 0. The results are obtained after extrapolations to the 
thermodynamic limit by finite-size scaling. 
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FIG. 2. Fourier transforms n(fc) of the charge density n c (x) 
for hole concentrations of n c h — 0.0625 and n c h — 0.125 in 
the conduction band and — Jh = Wt = 20J = 20J<j. Inset: 
charge density n c (x) for n% = 0.125. 

merical precision of the DMRG (\cr(i, j)\ < 10 _6 Vi, j). 

Next we investigate correlation functions by the 
DMRG method and start with the particle density in 
the conduction band n°{x), which shows Friedel oscil- 
lations induced by the open boundary conditions. For 
small hole densities n c h = 1/16, and ni = 1/8 we only 
find Ln c h j2 hole pockets (see inset Fig. indicating hole 
pairing. For larger hole densities however, a weak struc- 
ture in the hole pockets suggests that the two holes bound 
to a pair are separated by several lattice spacings. The 
pair binding energy in the low hole density region can 
be estimated by A pa ; r = 2Ei - E - E 2 « 0.016i for 
—Jh = lOt = 20J = 20Jd, where E n is the groundstate 
energy with n holes, and confirms pair formation. Also 
note the very weak amplitudes of the Friedel oscillations. 
In the Fourier-transform n(q) we only find one peak for 
very low hole doping ni — 1/16 at k = n^ir and two 
peaks at k = n^n and k = 2n^7r for larger hole dop- 
ing n c h = 1/8. In general, a Tomonaga-Luttinger liquid 
shows 2kp and Akp fluctuations, and 2kp — (n c + l)ir 
for a large Luttinger volume involving the electrons in 
the lower band. Our observation is consistent with these 
2kp and Akp periods for a large Luttinger volume. By 
fitting to the Friedel-oscillations for an impurity potential 
SrAx) oc C 1 cos{2k F x)x- { - 1+K "^ 2 + C 2 cos{4:kFx)x- 2K » 
p8[ we have determined the single correlation exponent 
K p w 0.51±0.05 for n° h = 1/8, indicating dominant CDW 
correlations. Because of trapped states near the bound- 
aries, these sites have to be discarded for the fit and the 
uncertainty stems from the fitting ambiguity. 

Independently, the correlation exponent K p can also be 
determined from the pairing correlations Pi{f)P} +r (f), 

where P/(/) = ^(c| jT c| +/4 =F cJ 4 c I+/,t) is either the 
singlet (— ) or triplet pair (+) creation operator. The 
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FIG. 3. Singlet pairing correlations Pi{f)P} +x {f) with 
20 J = 20 J d and n c h = 1/8. 

singlet pairing correlations decay as Pi(f)P} +x {f) oc 
A ln(.T) _1 - 5 x _1_1 / A "f + A 2 cos{2k F x)x^ K ^ 1/K " , and in 
agreement with the previous estimate we obtain K p w 
0.51 ± 0.05 by simultaneously fitting to / = 2,4,6, and 
8 (see Fig. ||). Note that the largest amplitudes of the 
pairing-correlation are obtained for / = 6 and 8, in agree- 
ment with the pair structure in n(x). 

Finally we turn to the spin-spin correlations shown 
in Fig. [|. As for the undoped case, all spin-spin cor- 
relation functions seem to decay exponentially SfSf +x oc 
cos(2fci?2;)e~ :E /'', and the correlation length £ = 7.9 ±0.1 
for n c h — 1/16 and £ = 11.2 ± 0.2 for n% = 1/8 remains 
close to the undoped case £ « 6.1. A more conclu- 
sive quantity is the string order parameter [ fL9[ g(x) — 

(IK. .,''- ) (E/=i,2^,/)) 

characterizing the Haldane S — 1 chain. The finite value 
of g{x) upon doping suggests that the spin- liquid state 
remains intact. However, for the doped case it is not 
clear whether g{x) remains finite in the thermodynamic 
limit, because our system sizes do not seem to be long 
enough if polaron-polaron correlations would destroy the 
string correlation. 

We finish by summarizing how all the findings includ- 
ing the surprising coexistence of gapless spin-modes with 
exponential spin-correlations (at intermediate distances) 
can be reconciled within the polaron picture. The holes 
which create a small fm cloud by the DE mechanism 
are bound to pairs as we see from the charge density, 
and 2kp and CDW correlations are dominant with 
K p w 0.51. The amplitudes of the CDW are very weak 
(see Fig. ||). This and the finite string-order suggest that 
the polarons are only a weak perturbation of the spin- 
liquid background of the 8 = 1 chain. Therefore we pro- 
pose a new model with two very different energy-scales 
for the spin-excitations. The larger one is of the order 
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FIG. 4. Spin-spin correlations SfS^ measured from the 
center and end of the chain for different values of hole dop- 
ing with — Jh = 10£ = 20J = 20Jd. Inset: string-correlation 
function g(x). 



of A s and stems from the spin-liquid background, giving 
rise to exponentially decaying spin correlations at the 
distances we can observe by the DMRG. The second en- 
ergy scale is very small and consists of the af interactions 
among the polaron-pairs. These are the gapless, lowest 
lying spin excitations, and at very long distances we ex- 
pect power-law decay of the spin-correlation functions. 
Since these interactions are among polaron-pairs, even 
our largest system-sizes contain only few such pairs and 
are too small to discriminate between exponential and 
power-law decay. Note however that £ S pin increases with 
doping. In order to estimate the energy scales we finally 
show the specific heat cv calculated by the TDMRG in 
Fig. ^. From the position of the large broad peak near 
T ~ 0.3i we see that the energy-scale of the spin-liquid 
remains almost unchanged upon doping, and from the 
sharp increase of cy at the lowest temperatures we esti- 
mate magnetic interactions among polarons at an energy 
scale of T ~ 0.02 J. 
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